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Abstract. We prove that the vorticity or the expansion vanishes for any shear- free 
perfect fluid solution of the Einstein field equations where the pressure satisfies a 
barotropic equation of state and the spatial divergence of the electric part of the Weyl 
tensor is zero. 
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1. Introduction 



There has been much effort devoted to establishing the shear-free fluid conjecture, that 
general relativistic, shear-free perfect fluids which obey a barotropic equation of state 
^ ■ p = p({i) such that p + fi 7^ 0, are either non-expanding (6 = 0) or non-rotating (u> = 0). 

This conjecture appears as such for the first time in Treciokas and EllispT] and has been 



established in many special cases, but a general proof or counter-example is still lacking. 
In support, the conjecture is known to hold in cases where: p = const (dust with a 
cosmological constant) [TT | f!2 | [18] : spatial homogeneity [T| [13]: dp/d/i = 1/3 (incoherent 
radiation) [2T] : dp/d/i = — 1/3 [HI [27] or 1/9 [26]; uj and ii are parallel [8]; vanishing of 
the magnetic part i?[Z] or of the electric part ^[HJ |9] of the Weyl tensor; 6 = 9(ji)fT5\ 
or 9 = 8(u)lQ\; Petrov types N|2j and III(3l H]; there exists a conformal Killing vector 
parallel to the fluid flow it [6]. It is noteworthy that there exist Newtonian perfect fluids 
with a barotropic equation of state, which are rotating, expanding but non-shearing. 
Hence, if true, such behaviour of fluids would be a purely relativistic effect. Recently, 
in an attempt to generalize the result established by Collins that H = u6 = 0, 
we managed to prove the conjecture for the case when divH = (that is when the the 
magnetic part of the Weyl tensor is solenoidal), with a 7-law equation of state[23]. This 
result was further generalised in a sense with our being able to establish the conjecture 
for the case when div H = and there is a sufficiently generic equation of state [5]. 
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In the present article we will focus our attention on generalising the result [TU E] that 
the conjecture holds when the electric part of the Weyl tensor vanishes, by considering 
spacetimes for which the electric part is solenoidal. Specifically, we prove 

Theorem Consider any shear-free perfect fluid solution of the Einstein field equations 
where the fluid pressure satisfies a barotropic equation of state and the spatial divergence 
of the electric part of the Weyl tensor is zero. Then either the fluid is non-rotating or 
non-expanding. 



2. Notations and conventions 



We shall be examining shear-free, perfect fluid solutions of the Einstein field equations 
Rab - ^Rg a b = H>u a u b + ph ab , (1) 

where u is the future-pointing (time-like) unit tangent vector to the flow, /x and p are 
the energy density and pressure of the fluid, respectively, and h ab = g ab + u a u b is the 
projection tensor into the rest space of the observers with 4-velocity u. The vanishing 
of the shear can be expressed by 

U a -b = -^Ohab + bJab - U a U b . (2) 

We shall assume familiarity with the notation and conventions of the orthonormal tetrad 
formalism as given by MacCallum [17]. We begin the analysis by choosing, what we 
believe to be, a well suited tetrad alignment for the physical problem at hand. This 
is the same as in our previous article [23]. First, e and e 3 are aligned with u and u>, 
respectively, such that uj = we 3 ^ 0. The relevant variables then become /i, p, 9, u, 
it a , Q a together with the quantities n a/3 and a a . Latin indices will be tetrad indices 
taking the values 1, 2, 3, 4. Greek indices take the values 1, 2, 3, while uppercase Latin 
indices the values 1, 2 and any expressions involving these have to be read modulo 3 or 
2 respectively. The sum of matter density and pressure will be written as 

£ = li + p. (3) 

Secondly, it is always possible, using the remaining rotational freedom and making use 
of the Jacobi identities and field equations, to further specialize the tetrad so as to 
achieve 

fij = Q 2 = Q 3 + uj = 0, (4) 
nu = n 22 = n. 

Herewith the tetrad is fixed up to rotations e.\ + ie 2 — > e %a {e.\ + ie-i) satisfying d§a = 0. 
Because of computational advantages, we will be replacing n a p (a ^ (3) and a a with the 
new variables q a and r a defined by n a _i Q+ i = (r a + q a )/2, and a a = (r Q — q a )/2. We 
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will also introduce extension variables, z a and j which are related to the components of 
the spatial gradient of the expansion, by 

d 9 = -9 2 /3 + 2u 2 + 3p)/2 + j, 

d a 9 = z a , (5) 

j being the (3+1) covariant divergence of the acceleration, 

J = u^ a = d a u a + ii a u a -u a (r a -q a ). (6) 

In order to relate the components of E ab with the spatial gradient of the acceleration, 
we will make use of the Ricci identity [T6] 

E(ab) = D {a U b) - LO( a U b ) + U( a U b ), (7) 

where S(ab) stands for the spatially projected and trace-free part of a tensor S a b- 

The complete set of initial equations of the formalism are now the Einstein field 
equations and the Jacobi equations, which we present, using the simplifications above, 
in the appendix. First notice that the equations ( 15711601) immediately lead to evolution 
equations for the variables r a and q a , 

3<V Q = -z a -9(u a + r a ), (8) 
3d q a = z a + 6(u a + q a ), (9) 

while ( 1501) and the (Oa) field equations ( I62ll64|) give us the spatial derivatives of u, 
2 

d t u = -z 2 -oj(q 1 + 2u 1 ), (10) 
2 

d 2 uj = -~zi + u(r 2 - 2u 2 ), (11) 

d 3 u = u(u 3 + r 3 - q 3 ), (12) 

together with the algebraic restriction 

2 

^33 = -Z-Z3. (13) 
3uj 

The evolution equation for n follows from (|58|) : 

Q 

d n = --n. (14) 

Acting with the commutators [d , d a ] and [di, d 2 ] on the pressure and using (1551) together 
with the conservation laws 

<9o/z = —89, d a p = —£u a , (15) 
leads to a first set of evolution equations for the acceleration and vorticity: 

d u a = p'z a — G9ii a , 

1 (16) 
d u = -u9(-2 + 3p'), 

where we have defined 

(17) 
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The spatial derivatives of the acceleration can be obtained from ((7j), using 

9aUa = ~\u 2 + - UA+iqA+l + UA-lfA-l - u\ + E AA (18) 

2 1 

<9 3 « 3 = g^ 2 + -j - iixqx + u 2 r 2 -uj + E 33 (19) 

diu 2 = -p'u6 + q 2 ui + ^n 33 ii 3 - «i« 2 + E 12 (20) 

<9 2 Mi = p'u;6> - nu 2 - ^n 33 u 3 - u x u 2 + E 12 (21) 

9iM 3 = ~u 2 n 33 - ^3^1 - M1M3 + E 13 (22) 

^2^3 = 7^i72 33 + g 3 w 2 - U2U3 + E 23 (23) 

9 3 «i = -^ii 2 n 33 + nu 2 + gi« 3 - + £ 13 (24) 

<9 3 w 2 = 7^1^33 - nui - r 2 w 3 - m 2 m 3 + £ 2 3- (25) 

Next we act with the [do, d a ) commutators on u and 6 and use the propagation of (fi~3~]) 
along u in order to obtain expressions for the evolution of z a along e$ and for the spatial 
gradient of j: 

d oZl = 0(-l + p>i - -w(-l + 9p')z 2 + -6cu(9G - 2)ii 2 (26) 

d z 2 = 0(-l + p')z 2 + + 9p')*i - -9u(9G - 2)ui (27) 

d z 3 = #(-!+ P>3, (28) 



ftj = - ^(27p' + 13)z 2 + ^i(18w 2 + 6 2 - 3j - 3fi) - —in 
o 3 2p' 

+ -#w(9G - 2)w 2 + 4w 2 gi (29) 

d 2j = p'6z 2 + \u(27p' + 13)zj + ^ 2 (18u; 2 + 2 - 3j - 3/i) - ^ 
6 3 2p' 

- ^w(9G - 2)ii x - Au 2 r 2 (30) 

dsj = p'6z 3 + \u 3 {6 2 - 18a; 2 - 3j - 3//) - ^U 3 - 4(r 3 - g 3 )u; 2 . (31) 
3 2p 

Now we may evaluate ^2 a [do, d a ]u a , using dSJ) and ( 1X81(19]) . which leads to an expression 
for the evolution of j in terms of d a z a . With the aid of the field equations fl68H70p . their 
propagation along u and the use of the [d , d a ] commutators on r a ,q a ,u a (excluding 
r i, 92> M3 ) together with the [di, d 2 ) commutators on uj and Jacobi equation fl53|) , 
one obtains algebraic expressions for the directional derivative d 3 z 3 , as well as evolution 
equations for j as 

_ . 6(G'£+p' -2Gp') . 2 M o „, . 

do] = — u + (1 - 2G)(z 1 u 1 + z 2 u 2 + z 3 u 3 ) 
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- ±9(l + 3G)j-p'9(l-9p')uj 2 (32) 
9(G'£-2Gp'), 2 2 , 2G . 

3Zs = P '(i + 3 P ') ^ ~ ^ ~ IT%r lMl + Z2U2 ~~ ^ ~~ 3 

2G + 9p' 2 -6p' + l g(- 2 + 3G) F 

+ : , Q , W + — — £ 3 3 + ^2 - 33 

1 + 3p' 1 + 3p' 

(G" = and -u 2 = wf + 'U 2 , + «§)■ Also the solutions of equations ( 153|) and (|65|) are given 
by 



9{2Gp'-G'£) f . 2 2G , . . . , 

^ = 3^(1 + W) {u ~ 3%) + Mi + W) (ziUl Z2U2 ~ 3Ms) 



and 



r 3 -3q 3 + 3u 3 9(-2 + 3G)„ 9p' 2 + 2G - 9p' 

H -23 ; ; ; 9uj 

3u 3u{l + 3p') 3(1 + 3p') 

- -^-(r 2 z 2 - qizi) + (n + qi)r 2 + n(q 3 + r 3 ) - E 12 (34) 

9{2Gp' -G'£) ,. 2 . 2 ~, 2G . . . n . , 

^ 2 = 3^(1 TW) K + " 2 " %) + 3c(l + 3 P ( ^ 1 + ^ " 2Z3%) 

3r 3 - g 3 + 3m 3 0(-2 + 3G) _ 9p' 2 + 2G - 9p' 

H z 3 r-^33 } ; 9bj 

3u 3 3u{l + 3p') Si 3(1 + 3p') 

- -^~(r 2 z 2 - qxzi) - (r 2 + q 2 )qi - n(q 3 + r 3 ) + #12- (35) 

From the Ricci identity [IB] . i?( a b) = 2u^ a c0f ) ) +D( a Ub), the components of the (trace- free) 
magnetic part of the Weyl curvature are given by 

H n = - u(ii 3 + r 3 ), #22 = -u(u 3 - q 3 ) 

Hi 3 = z 2 /3 - uqi, H 23 = -Zi/3 + ur 2 , H X2 = 0. (36) 

At this stage we introduce the condition that E is solenoidal, divi£ = 0, which by the 
Bianchi identity [16], (divi£) a = \D a n — 3u b H a b reduces to 

£ iii + 3p'u(z 2 — 3uqi) = 0, (37) 

£u 2 — 3p'co(zi — 3ur 2 ) = 0, (38) 

{£ + 18wV)«3 - 9wV(-r 3 + 93) = 0. (39) 

Apart from /icu9£ ^ we shall also assume that p' ^ 0, ±1/3, 1/9, u\ + u 2 ^ 0, E ab ^ 0, 
H a b 7^ 0, 9 not just a function of either fi or u, so as to avoid duplication of known 
results. Propagating equations fl37|) and fl38l gives then 

3p'[2£ - 9(1 - 3p')u 2 }z A - 9[2£(6p' - 2) + 81Gp'co 2 ]u A = 0, (40) 

These two equations show that z A is parallel to ii A , unless 

2£ - 9(1 - 3p')u 2 = 2£(6p' - 2) + 81Gp'cu 2 = 0. (41) 
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Taking the evolution of ii2Z\ — U\Z% = along u, this results in 

27p'(3G + 3p' - l)co 2 - £(-1 + 27Gp + 9p' - 54p' 2 ) = 0. 

If the coefficient of cu 2 were zero, the latter equation would give 
£(-l + 9p') 2 = 0. 

Hence, the vorticity is a function of the matter density and one has <9i/i<9 2 u;— <9 2 o;<9iyU = 0, 
which reduces to (iiiZi + U2Z-i)£ = 0. It follows that u\ = w 2 = 0, which is impossible. 

4. iiA not parallel to za 

Eliminating u from ( |4T|) fixes the equation of state to be 

9Gp' - 18p' 2 + 9p' - 1 = 0. (42) 

Note that a 7- law equation of state is then only possible when p' = —1/3 or 1/9, so 
henceforth we will exclude all sub-cases in which p' is constant. By (f4Tj) the vorticity is 
again a function of matter density. It follows that 

d A ^d 3 uj - d A ud 3 n = 0, 

which after eliminating g 1; r 2 and q 3 using equations fl37|) . f )38|) and fl39|) yield 

£(-3u;ua + za+i)m 3 = 0. 

It follows that u 3 is zero, as otherwise, propagation of equation f j4T1) along would 
imply £(— 1 + 3p') , ii J 4 = and again -u^ = 0. 

When ii 3 is zero, the expressions (HE]), $22]), ([23]), (EE]) and ([T3j) show that 
Z3 = 0, n 33 = 0, E 13 = «ir 3 , £ 23 = -w 2 g 3 , 
E33 = -g (2w 2 + j) + wigi - m 2 ^2- 



Moreover, the conditions div.E = 0, (I37H39P can now be rewritten as follows 
Sill + 3ojp'(z2 — 3wgi) = 0, 

£u 2 — 3up'(zi — 3ur 2 ) = 0, (43) 
P'(q3 - ^ 3 )w 2 = 0. 

Propagating equation (T4T!) (a) along ei and e 2 , using equations ( 143]) . yields 

£{-u x + 3A0 = 0, (44) 
£[ii2 + 3pV 2 ) = 0. (45) 

Again propagating the latter equations along e\ and e 2 respectively and eliminating E 12 
shows then 

99p' 2 (3p' - l)(18p' 2 - Ihp' + l)j + 9p' 2 (3p' - l)(36p' 2 - 9p' + 5)(«i2i + u 2 ^ 2 ) 
-90(3// - l)(162p' 4 + 108p' 3 - 162p' 2 + 21p' - + w 2 ,) 

+6p' 3 9(27p' 2 - 18p' + 7)£ = 0. (46) 
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Now we use z 3 = and equations (14"TH4"1) . (|4"5|) and (14"31) to simplify equation (1551 ) to 

(3p' - l)(162j/ 4 + 108p' 3 - 162p' 2 + 21p' - l){u\ + u\) - 6(27p' 2 - 18p' + 7)£p' 3 

-9p' 2 (3p' - l)(18p' 2 - IV + l)j = 0. (47) 

This can now be combined with ( 1461) to give 

p ,2 (3p' - l)(36j9 /2 - V + 5)(«iui + z 2 « 2 ) = 0. 

It follows that m,4 is orthogonal to 2^, after which propagation of equation (j46H along 
e gives 

2(1296p' 4 - 2106p' 3 + 522p' 2 + 9p' - 5){3p' - l) 2 (u\ + 

-2Ap' 2 8(p + l)(81p' 3 - 207p' 2 + 93p' - 7) = 0. (48) 

Taking the time derivative of the latter equation, eliminating j, z\ using equation (|46l) 
and z 2 U2 + ^1^1 = 0, results in a new equation between p,p' and m 2 + u 2 ,. Eliminating 
the acceleration between the latter pair eventually results in 

£p' 3 (3p' - l)(18p' 2 - I5p' + 1) 
1296p' 4 - 2106p' 3 + 522p' 2 + 9p' - 5 X 
(26244p /6 - 59778p' 5 - 106515p' 4 + 156249p' 3 - 43065p' 2 + 3267p' - 26) = 0, 

showing that p' has to be a constant, in contradiction with |42] and our assumptions 
p^ -1/3, 1/9. 

5. Discussion 

Newtonian homogeneous cosmologies are known [T9] to provide counter-examples to the 
Newtonian shear-free fluid conjecture. In these models the spatial divergence of the tidal 
field (the Newtonian analogue of E) automatically vanishes, such that, when searching 
for possible counter-examples to the relativistic shear-free fluid conjecture, it is tempting 
to look for candidates within the class of perfect fluids having a solenoidal electric 
part of the Weyl tensor. The result obtained in this paper shows that this attempt 
remarkably (or luckily — depending on one's attitudes towards the conjecture) fails: 
shear-free perfect fluids obeying a barotropic equation of state and having divi? = are 
all non-rotating or non-expanding, in contrast to their Newtonian brethren. Whether 
the conjecture in its full generality is true or false still remains an intriguing issue. 
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6. Appendix A 

Commutator relations, using a a p = 0: 

[do, <9i] = ii^o - 01&! + (w 3 + £l 3 )d 2 - (u 2 + ^ 2 )<9 3 
[do, d 2 \ = u 2 d - (u 3 + fi 3 )<9i - S 2 d 2 + (u l + fii)<9 3 
[d , d 3 ] = u 3 d + (u 2 + Q 2 )9i - 9 3 d 3 - (wi + ^i)d 2 
[<?i, <9 2 ] = -2w 3 (9o + g 2 <9i + ri<9 2 + n 33 <9 3 
[<9 2 , <9 3 ] = -2wi9 + g 3 <9 2 + r 2 d 3 + n n <9i 
[<9 3 , di] = -2u 2 d + r 3 di + q x d 3 + n 22 <9 2 

Using in addition the simplifications Ui = u 2 = Q x = Q 2 = 0, Q 3 = —u 3 = —ui, 
nu = n 22 = n one obtains the: 

Jacobi equations: 

d 3 u = uj(ii 3 + r 3 - q 3 ) (50) 

din + <9 2 r 3 + <9 3 g 2 - n(n - gi) - r 3 r 2 + g 3 g 2 = (51) 

<9 2 n + d 3 ri + <9ig 3 - n(r 2 - nq 2 ) - r 3 ri + q 3 qi = (52) 

2 

dir 2 + <9 2 gi + <9 3 n 33 - -9u - r 2 r\ + q 2 q x - n 33 (r 3 - g 3 ) = (53) 

d 2 u 3 - d 3 ii 2 - rail - <73«2 - r 2 ii 3 = (54) 

d 3 iii - diii 3 - nu 2 - r 3 iii - q x u 3 = (55) 

4 

2d u) + diii 2 - d 2 iii - iiiq 2 - ii 2 r 1 - u 3 n 33 + -9co = (56) 

o 

1 9 

d (r a - q a ) - -d a 9 + z a + -(r a -q a + 2u a ) = (57) 
3d n + 3d 3 cu + n9 — 3u(u 3 + r 3 — q 3 ) = (58) 
3d n 33 + 3d 3 u + n 33 9 - 3u(ii 3 + r 3 - q 3 ) = (59) 

d {r a + q a ) + U(r a + q a ) = (60) 
Einstein equations: 

d 9 = - l -9 2 + 2uj 2 - 1 -{^ + 3p) + j (61) 

2/3^1 + <9 2 w-u;(r 2 -2M 2 ) = (62) 

2/32 2 - diu - u{qi + 2ui) = (63) 

2/3z 3 - wn 33 = (64) 
-<9ir 2 + d 2 qi - r x r 2 - q x q 2 - 2r 2 q x - 2r 3 n + n 33 (r 3 + q 3 ) - 2g 3 n = 

-diu 2 - d 2 iii - 2uxu 2 + inq 2 - u 2 r\ (65) 
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-<9 2 r 3 + d 3 q 2 + d x n - din 33 - n(ri - <?i) - 2q 1 n 33 - r 2 r 3 - q 2 q 3 - 2r 3 q 2 = 

-d 2 u 3 - d 3 ii 2 + ui{n 33 - n) - 2u 2 u 3 + u 2 q 3 - u 3 r 2 (66) 

diq 3 - d 3 r 1 + d 2 n 33 - d 2 n - r 2 (2n 33 - n) - nq 2 - rir 3 - q x q 3 - 2r\q 3 = 

-diii 3 - d 3 iii - iii(2u 3 + r 3 ) + u 2 {n - n 33 ) + ii 3 qi (67) 

-8^ + + <9 2 g 2 - <9 3 r 3 + ?2 + r 3 + ^ n 33 ~ nn 33 + »"i + 5? ~ ^2^2 - r 3 (? 3 = 

2 8 1 

-6» 2 + -w 2 - -(2[i - j) - 9iui - m 2 + w 3 r 3 - w 2 g 2 (68) 

-<9 2 r 2 + 9 2 g 2 - din + 5 3 g 3 + g 3 + r\ + in 33 - nn 33 +rj + q 2 - r x q x - r 3 q 3 = 
2 8 1 

-O 2 + -w 2 - -(2/i - j) - <9 2 m 2 - m 2 + uxn - w 3 g 3 (69) 

-d 3 r 3 + <9 3 g 3 + diq x - 9 2 r 2 + q\ + r 2 - ^n 2 3 + r 2 + g 2 - r x q\ - r 2 q 2 = 
2 1 

-6 2 + -(2w 2 - 2/i + j) - 9 3 m 3 - ii 2 3 + - uigi (70) 
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